Cyclostationary measurement of low-frequency odd moments of current fluctuations 
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Measurement of odd moments of current fluctuations is difficult due to strict requirements for 
band-pass filtering. We propose how these requirements can be overcome using cyclostationary 
driving of the measured signal and indicate how the measurement accuracy can be tested through 
the phase dependence of the moments of the fluctuations. We consider two schemes, the mixing 
scheme and the statistics scheme, where the current statistics can be accessed. We also address 
the limitations of the schemes, due to excess noise and due to the effects of the environment, and, 
finally, discuss the required measurement times for typical setups. 



I. INTRODUCTION 

The discrete nature of the charge carriers shows up 
in the peculiar statistics of the transmitted current 
through mesoscopic systemsi Consequently, these statis- 
tics in general cannot be described through the usual 
Gaussian probability density which is completely deter- 
mined by its first two moments (cumulants), average 
current and noise, whereas its higher cumulants van- 
ish (for a difference between the moments and cumu- 
lants, see Ref. Recently, there have been many 
theoretical predictions for the behavior of either the 
third moment of fluctuations or the full couting statis- 
tics in different types of mesoscopic systems (see a review 
in Ref. 13), including normal-metal tunnel junctions^ 
normal-metal-superconductori^ and superconductor- 
superconductor junctions, ^'^'^ diffusiveiSiiiiiSiii and 
chaotic wires, double-barrier junctions^^ and Coulomb- 
blockaded systemsii^ For example, the second and 
third moments of current fluctuations through a phase- 
coherent scatterer at T = may be expressed as 



S = (SISI) 



R = (SISISI) 



e'\V\ 



^T„(l-T„) = F2e|J| 



^r„(l-T„)(l-2r„) = F3e2/, 

(1) 



where T„ are transmission eigenvalues. Their distribu- 
tion depends on the properties of the scatterer^ and 
it can be characterized by defining the Fano factors F2 
and F3. Compared to the average current J, the higher 
moments may thus reveal additional information about 
the scatterer (for an example, see Ref. [T^ . Moreover, 
R implies that the distribution of fluctuating cur- 
rents is " skew" . This has important consequences in the 
situations where the driven current fluctuations act as an 
environment to another mesoscopic system^ 

However, up to date there exists only one published 
measurement of the higher (than second) moments^^*^ 
One of the reasons for this is the necessary conditions 
for the filtering in these devices: fmpedance matching 
between the amplifier and the sample is possible typi- 
cally only within a fairly narrow frequency band, outside 



of which the signal from the sample does not couple to 
the measurement device. However, as we show below, 
for the typical measurements of the odd n'th moments 
of fluctuations, the requirement for the bandwidth 26lu 
of the measured frequencies around the mean frequency 
ujo is Suj > LUo/n^ Especially for the measurement of 
high-impedance samples (Z larger than 1 kfl) this re- 
quirement is very hard if not impossible. In this paper, 
we show how this requirement can be circumvented by 
a cyclostationary driving of the measured system (see 
Fig.lIJ. 
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FIG. 1: Measurement schemes considered in this paper. Each 
contain an input / filtered through a band-pass filter (denoted 
with BP, corresponding to the function Hbp). In the "DC 
mixing scheme", analogous to that utilized in Ref.l20l the sta- 
tionarily driven and band-pass filtered input is mixed twice 
with itself, low-pass flltered (denoted with LP, corresponding 
to the function H^p) and averaged. The "AC mixing scheme" 
is a modification of the DC scheme, including a cyclostation- 
ary driving of the measured system, and mixing with an extra 
function f{t) — cos(na;o(i + to)), synchronized with the driv- 
ing and with a controlled phase shift to- In the "AC statis- 
tics scheme" , the statistics of the output is measured and the 
driven and flltered signal is mixed only with f{t). 

To motivate the use of cyclostationary driving, con- 
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sider the n'th central moment of current fluctuations in 
the frequency space 

M(")(wi,C^2,...,^„) = {5I{uJl)SIiLU2)...6I{LOn)), (2) 

where the brackets denote ensemble averaging and 
61 (uj) = I {to) — (/(w)). If the studied system is time 
independent, i.e., if it is driven with a constant voltage, 
we get the signal only from the frequencies that satisfy 

LJl + UJ2 + ■ ■ - ^n ^ 0- (3) 

Consider now a situation where the signal is band-pass 
filtered before its correlators are measured, such that only 
the frequencies in the band uj G [wq — + Suj] pass 

through the filters. If 6oj ^ ujq, we cannot fulfill con- 
dition (PJ for odd n. Simple algebra shows that for the 
n'th moment (where n is odd), we only get some signal 
if 

Suj > oJo/iT'. (4) 

However, this condition is true for stationary signal, but 
if we adiabatically drive the system with frequency lod, 
condition Q changes to 

U>[ + LUi + . . . Un-l — mUJD, (5) 

with an integer ni. Reasonable examples considering the 
odd n'th moments are lod = ij^o, Swq, . . . , nujQ. 

In the second section of this paper, we detail two spe- 
cific measurement schemes using cyclostationary driving 
and compare them to that applied in Ref . "^o". The first of 
the schemes is based on mixing the signal with itself and 
with an oscillatory function f{t), and the second on mea- 
suring the statistics of the output signal from the sample 
after filtering and mixing with f{t). The third section 
discusses the effect of cyclostationary driving in some ex- 
ample systems and the fourth considers the limitations 
of the presented measurement schemes due to the excess 
noise in the environment, and estimates the required av- 
eraging time for obtaining an accurate signal. The details 
of the calculations are presented in the Appendix. 



II. MEASUREMENT SCHEMES 

Let us consider the specific measurement schemes for 
measuring the moments of current fiuctuations depicted 
in Fig. The emphasis of this analysis is on the 
third moment, but the results and the schemes are fairly 
straightforward to generalize also to higher moments. In 
both of these schemes, the studied system is driven with 
a time-dependent voltage of the form (for most of the 
discussion below, the different harmonics do not need to 
be in the same phase, but we assume them in phase for 
simplicity) 

V{t) =Vo + Vi cos(woO + cos(3woi)- (6) 



This voltage produces in general an output current of the 
form 

I{t) =^Iicos{lujot). (7) 
I 

In linear systems with resistance R, Ii —Vi/R and thus 
Ii = for / > 3. In nonlinear systems, the output current 
would also contain higher harmonics of wq, but these do 
not contribute essentially to the results. However, the 
coefficients Ii should be taken from the proper Fourier 
analysis of the output current. 

With pure voltage driving, the moments of the low- 
frequency current ffuctuations are proportional to either 
the average current I{t) (for odd moments) or its absolute 
value (for even moments). This holds as long as the driv- 
ing is adiabatic, i.e., the highest driving frequency 3ujq is 
low compared to the internal energy scales of the probed 
system or to eV{t)/h,^^^ In fact, any slow measurement 
of the signal as a function of the driving voltage can be 
considered as "adiabatic" in this sense — the only dif- 
ference to that considered here is that typically the time 
scales of such voltage variations are of the order of sec- 
onds, much slower than the scales for band-pass filtering. 
Furthermore, for low frequencies, the fluctuations can be 
considered "white", i.e., they are of the form (below, we 
reserve the symbol S for the second moment, i.e., noise, 
and use the symbol R for the third moment )2S 

S{t, t') = {5I{t)5I{t')) = ^Si cos{luoQt)5{t ~ t') (8) 

I 

and 

R{t,t',t") = (SI{t)SI{t')SI{t")) 

= ^ i?i cos{lLJot)S{t - t')d{t - t"). (9) 
I 

Similarly, the n'th central moment would be 

M(")(ii,t2, . . . ,t„) = {dI{ti)SI{t2) . . . SI{tn)) 

^ J2 ^^/"^ cos{lu;oti)6iti - t2)5{ti -ts)... 6{ti - i„)- 

(10) 

In most cases below, it is enough to cut the series in the 
n'th harmonic (for the n'th moment) as the higher ones 
do not contribute to the results. There are many sources 
for the higher harmonics than those present in the driv- 
ing and one has to always first do the proper Fourier 
analysis to find these, prior to examining the results: i) 
finite temperature of the order of the voltage V makes the 
relations S{V) and R{V) nonlinear even in many linear 
systemSfiiS ii) if the sign of the total current is allowed 
to vary, one has to take the Fourier components from the 
absolute value of the current for S{V), and ni) these har- 
monics arise naturally in nonlinear systems. Examples of 
these are discussed in Sect. IIIl| However, we stress that 
all these effects can only make rise to frequency compo- 
nents of the form Iujq^ where I is some integer. 
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For what follows, the assumption of the "white" fluc- 
tuations is in fact less strict than assumed in Eqs. (|8I10|I : 
it is enough that the fluctuations are frequency indepen- 
dent in the window allowed by the band-pass filter. 

For both schemes, the signal is assumed to be filtered 
with an ideal band-pass filter allowing frequencies u; S 
[loq — 5uj^ujq + Slo] . Such a filter can be described through 
the boxcar function Hbp{'^) = 0{\uj\ — {luq — Suj)) — 9{\uj\ — 
{ujQ + Suj)), where 9{uj) is the Heaviside step function. 
Thus, the outgoing signal from the filter is in the time 
domain a convolution (below, unless it is obvious, the 
functions in the frequency domain are denoted with a 
tilde) 

X{t)^^HBPit)^m. (11) 

V ZTT 

At the end of the processing, the signal is gathered for a 
time exceeding I/ujq. Such a process can be described by 



a convolution with a low-pass filter H^p{uj) = l — 9{\oj\~ 
LUjn), where aj^n *C ^o- 



A. Mixing schemes 

In the mixing schemes depicted in Fig. ^ (schemes a 
and b), we assume that the signal is, after filtering, twice 
mixed with itself and with a function f{t) specified below, 
low-pass filtered, and finally ensemble averaged. Thus, 
the measurement result is 

{Y)^^{HTp^[fit)XHt)]). (12) 

Assuming that over a sequence of measurements per- 
formed for the averaging, the filters and the function f{t) 
do not essentially change, one obtains 



47r' 



-H'rp,[f{t'mBP^m] 



fit') I dhdt2dhHBp{t' ~ ti)HBp{t' - t2)HBp{t' ~ t3){I{ti)I{t2)I{h)) 

I 



(13) 



Below, we take a closer look at two ways of driving the 
signal: first with a DC bias and wide band-pass filters, 
and then with an AC bias and narrow band-pass filters. 



1. DC mixing scheme 

Assume first that the system is DC biased, i.e., only 
the first term in Eqs. H6I9I) is nonzero. In this case, we 
do not need the function generator, i.e., f{t) = 1. Such 
a setup corresponds to the one utilized in Ref . and it 



is depicted in Fig. 

One can write the current in the form I{t) = {I)+SI(t), 
where 51{t) is a zero-mean fluctuating current. Then the 
third raw moment of current fluctuations is 

{I)S{S{t - t') + d{t - t") + S{t' - t")) + R5{t - t')6{t - t"). 

(14) 

The average signal coming out of the sampling filter H^p 
is 



{Y)^—{Hl^P^{HBP*I{t)f) 



^HTp * I dt'idt'^dt'iHBPit' t'i)HBp(t' t'i)HBp{t' tl){I(t'i)I(t'i)I(t'i)) 
(/) j dt"HBp{t' -t")j + i{I)sl^j dt"HBp{t' -t")^ l^j 



dt"Hlp{t' -t"] 



R j dt"Hlp{t' - t") 



(15) 



I 

The band-pass filter filters out the DC signal, as long as it is a band-pass filter, i.e., 5u) < loq. Therefore, 

/ dt"HBp{t - t") V2^Hbp{0) = (16) 
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only the last term remains, and we obtain 

{Y)^R^{i5ijj-LOQfe{i5uj-LOo). (17) 

It can be seen that for 5uo < wq/S, {Y) = 0, in agreement 
with the discussion in Sect.Q] The scaling of the observed 
quantity as a function of the bandwidth was applied in 
Ref. to show that it was indeed the third moment that 
was measured. 



2. AC mixing scheme 



Assume that the current, noise and the third moment 
are of the form given in Eqs. Q, and ©, respec- 
tively. For simplicity, we assume that the lowest driving 
frequency equals the center frequency of the band-pass 
filter and denote both by ojq. The average output signal 
is given by 



(Y) -i^rp * \ fit') 



j dt"HBp{t' - t") cos(/wot") j 
^(^IiJ dt"HBp{t' - t") cosilujot")^ Si J dt"Hlp{t' ~ t") cos(/wo<") j 
j dt"Hlp{t' - t") cos{lujQt") 



(18) 



Now we have three types of integrals with the band-pass filter Hbp ■ We obtain 
dt"HBp{t' - t") cos{lLJot") = V2TTHBp{liJo) cos{lujot') = V27r5i,i cos(wot') 
dt"Hlp{t' - t") cos{lLUQt") = {Hbp * Hbp)\uj=Iujo cos(/woi') = 2(5w [2(5/,o + Sia cos(2woi')] 



(19) 
(20) 



dt"Hlp{t' -t")cos{lujot") 



Combining these results, we get 
H 



{Hbp * Hbp * Hbp)\uj=Iujo cos(Zwot') 



'LP 



fit') 



27r/f cos^(wot') -I- GhSuj [2Sq cos(woi') + 5*2 cos(2woi')] 



v/2^ 



2tt 



[3(5;4 cos(woi') + <^;,3 cos(3woi')] 



(21) 



[3i?i cos(ti;oi') + R3 cos(3ti;o^')] 



r 



Now it is time to specify the mixing function f{t). This is dependent only on the third moment R3. Moreover, 
mixes some of the above oscillating functions to the the dependence on the phase shift shows that the func- 
zero frequency and thereafter the low-pass filter i?™p tion f{t) and the driving have to be synchronized. This 
filters the other frequencies out. There are two mean- phase dependence can also be utilized as a signature of 
ingful choices for f{t), (a) f{t) = cos(ti;o(^ + ^o)) and the fact that the measured quantity really corresponds to 
(b) f{t) = cos{3L^o{t + to)), where to represents the phase the properties of the signal coming from the sample, anal- 
shift between the function generator and the signal. With ogously to the dependence on the bandwidths utilized in 



the first choice, the result is 



Ref. 



(Ya) 



Ri 



cos(woio), (23) 



AC statistics scheme 



and with the second, 
'1 



(n) = ( + ^i?3 ) cos(3c^oio). 



(24) 



In the second scheme (Fig. ^) , the AC-driven current 
is mixed only once with the function f{t) — cos(wo {t+to)) 
and the output signal 



This result illustrates that, in a linear system, setting 
Iq = Ii = h = 0, the signal with f{t) = cos(3wo(t -I- ^o)) 



Y - —H^^p * [/{Hbp * /)] 



(25) 
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is measured many times. This produces an ensemble of 
y-values, whose moments are related to the moments of 
the current I{t). As above, we assume that the band- 
widths of both the band-pass filter, Slo, and of the sam- 



pling, iVrn, are much lower than loq. 

The n'th central moment of this y-distribution is given 

by 



lldUHrp{U)f{U) 



.i=l 



.4=1 



(26) 



where I labels harmonics of the base frequency luq. Utilizing the white- noise assumption as in Eqs. (|8I10|I yields 



l[dUHrp{U)f iU) 



dt' 



\{HBp{U-t') 



m/"^ cos{lu;ot'), 



(27) 



(n) 

where is the coefficient of the Pth harmonic of the n'th central moment in the measured signal. This integral is 

evaluated in Appendix^ The results for the six lowest moments are (the first, average, is the raw moment, the rest 
are central moments, i.e., defined w.r.t. the average) 



(Y) = cos{ujoto) 



iiSYf) 
{(SYf) 



An 
327r2 

, ,3 



[25o + S'2COs(2u;oio)] 

[3i?i cos(cjoio) + -R3 cos(3a;oio)] 



247r3 
61447r4 
12807r5 



QM^^^ + AM!f^ cos{2LUoto) + Aff cos(4woio) 



lOMf ^ cos{ujQto) + 5A/f ^ cos(3cjoio) + M^'^^ cos{5ujQto] 



20A'4®' -I- IbA^"' cos(2woio) + 6Mf ^ cos(4woio) + M^"' cos(6a;oio) 

I 



r(6), 



r(6). 



r(6), 



(28) 
(29) 

(30) 
(31) 
(32) 
(33) 



where LUd — min(wm, Slo). Note that which harmonics are 
nonzero depends on how the system is driven. For exam- 
ple, if one drives the system with a voltage V3 cos(3ti;o0i 
the only finite moments (among the first six) contribut- 
ing to the measured signal are 5*0, i?3, Mq'^\ Mq^^ 

and Mg^' . Such examples are discussed in more detail in 
the following section. 



III. EXAMPLES OF HIGHER HARMONICS 



In this section, we discuss examples of systems where 
one could measure the third moment of fluctuations 
with AC driving. The considered examples are the 
simplest linear systems at low temperatures, the finite- 
temperature case where part of the signal is mixed to 
higher harmonics, and thirdly a generic example of a non- 
linear system. 



A. Linear systems at low temperatures 

The simplest but yet important example of a sys- 
tem where our theory is applicable is a linear system 
with resistance R, driven with an AC voltage V{t) — 
I^Rcos{3LUQt). In a tunnel junction, neglecting the ef- 
fects of its environment f22i2^ the third moment is inde- 
pendent of the temperature^ whereas the noise is in gen- 
eral a mixture of shot noise (produced by driving) and 
thermal noised In other systems, also the third moment 
becomes temperature dependentii^ For eV kpT, the 
output current, noise and the third moment are 



I{t) — I3 cos{3ujot) 

S{t) = eF2h\cos{3ujot)\ = 6^2/3 (- + ^ cos(6cjoi) 

TT OTT 



R{t) = e^i^3/3Cos(3tJoO- 



(34) 

(35) 
(36) 



Here F2 and F3 are the Fano factors {F2 = F3 = 1 for 
a tunnel junction) for the second and third moments, 
respectively. Applying the AC mixing scheme and mixing 



6 



with f{t) = cos(3wo(i + ^o)) yields the output 

In the AC statistics scheme, the average of the output 
signal vanishes, the variance is 

{{oY) } = 6^2/3 (38) 

and the third moment is 

iiSY)^) = I (^^^^%iM)'cos(3c.oto)e^F3J3. (39) 

This illustrates that, as long as the thermal noise and 
environmental noise are negligible, the best signal from 
the third moment (compared to the second) is obtained 
with wide-band filtering but low-amplitude driving. 

For a specific example, assume a realistic set of num- 
bers, h = 100 nA and min(w„i, (5w)/(27r) = 10 MHz. 
Then, we get 

~ i^2(320pA)2 (40) 
« i^3(46pA)3 cos(3a;oio) (41) 

Thus, the third moment should be well observable in this 

situation. With larger currents or smaller bandwidths, 

the second moment becomes relatively larger, but only 
as ((y2))i/2/((y3^)i/3 ^ (/3/(e5^))i/6. 

B. Diffusive wire at a finite temperature 



Here Zs is the resistance of the measured sample. 

For example, let us assume that we drive the system 
with the voltage 

V{t) ^V^cosi^LUot). (48) 

As a diffusive wire is a linear system, the average current 
follows the oscillations of the voltage, I{t) = I3 cos(3a;ot). 
Our measuring schemes yield information about the 
stationary part Sq of the noise, and about the Sluq- 
component i?3 of the third moment. These are plotted in 
Fig- Has a function of eV3/2kBT. One can observe that 
in this particular case, the crossover to the "pure shot 
noise" takes place for eVs ^ lO/csT whereas the crossover 
between the low- and high-temperature third moments is 
wider, saturating only at some eVa AOksT. 

0.45 1 • • • • • 1 

0.4 '; 

30.35 




5 10 15 20 25 30 

eV3/(2k3T) 



The temperature dependence of the higher moments of 
fluctuations depends on the studied system. For example, 
for a diffusive wire one expectsi^ 

S = — coth + —2— (42) 

3 \2kBTj 3Zs ^ ' 

6(-l + e4p) + (l-26e^P + e4p)p 
R = el ; TT—^ • 43 

Here p — eV/2kBT. In this case, for an oscillating volt- 
age, some of the moments mix into higher frequencies, 
altering the expected outcome of the measurement. How- 
ever, in the extremal cases eVi <C ksT and eVi 3> ksT, 
where eVi is the amplitude of the oscillating voltage, this 
mixing is not important and one may simply use 



S{e^ < ksT) 



2kBT 



S{eV, » kgT) - -e|/,(t)| 

R{eV, « ksT) ieV,(i) 

RieV^ » fcsT) = ^e^V^it). 

15 



(44) 
(45) 
(46) 
(47) 



FIG. 2: The stationary part So - 4fcflr/(3Zs) of the sec- 
ond moment (dashed line) and the lowest harmonic Rs of 
the third moment (solid line) in a diffusive wire driven with 
V3 cos(3a;ot), normalized to e/3 and e^Is, respectively. Note 
that the previous diverges with eVs <^ ksT as then So con- 
tains only the thermal noise and is essentially independent of 
/3. 



C. Generic nonlinear system 

Many interesting mesoscopic systems are inherently 
nonlinear: there exists an energy scale Ec below and 
above which the electron transport mechanism may dif- 
fer. Such a behavior is often signalled by different Fano 
factors for the second and third moments below and 
above Ec- Below, we explain how the moments can 
be measured in such systems applying cyclostationary 
driving. There are two extreme alternatives: either one 
drives the system with a single oscillating voltage, say, 
V{t) — V3 cos(3a;oi), or with a large DC voltage Vb fixing 
the operating point and with a small AC excitation of 
the form V3 cos(3ci;oi) on top of it. The first alternative 
has the benefit of minimizing the excess noise compared 
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to the third moment (important for the averaging time), 
whereas in the latter case the results are easier to inter- 
pret. 

In nonlinear systems, the nonlincarities mix the driving 
frequencies to higher harmonics of the base frequency. To 
see an example of this, we consider a generic non-linear 
system at zero temperature, assuming that the differen- 
tial conductance of the system is 



dl 
dV 



Ga, V <V, 
Gb, V > Vc. 



(49) 



Here Vc is some characteristic voltage describing the non- 
linearity. For example, in the case of Coulomb blockade, 
it would be determined by the charging energy, or in 
the case of a normal-metal - insulator - superconductor 
(NIS) junction, by the superconducting gap. Typically in 
such systems, the process for charge transport below and 
above Vc differs (for example, in the case of intermediate- 
transparency NIS junction, for V < A/e, it would be due 
to Andreev reflection, and for V > A/e, due to quasipar- 
ticles). Therefore, also the differential Fano factor may 
differiSS. between the two regimes. Thus, let us assume 
that the differential noise (second moment) is given by 



dV 



eF^GaSgniV), V < Vc 
eF^GtSgn{V), V > Vc- 



Similarly, the third moment would be 



V <Vr 




e'FiGb, V > Vc. 



(50) 



(51) 



Given a slowly varying voltage V{t), the resulting second 
and third moments would be 



Sit) 



R{t) = 



eF^^Ga\V{t)l \V{t)\ < Vc 

e\F^GaVc + F^GbiVit) - K)|, \Vit)\ > Vc 

(52) 

'e^F^GaVit), \Vit)\ < Vc 
e^F^GaVc + F^Gb{V{t) ~ Vc), \V{t)\ > Vc. 

(53) 



Now drive the system with V{t) = V3 cos{3ujQt), V3 > Vc. 
The DC part of the noise, contributing to the measure- 
ment outcome, is 



2eV. 



FSGa + {F^Gb~F-Ga)\ll~[ ^ 



(54) 

The lowest harmonic of the third moment, oscillating 
with the frequency SujQt would then be 



i?3 = e^F^GaVs 



X V3 arccos 



2e^ 

TT 

Vc 
V. 



{F^Gb - F^Ga) 



Vc\l- 




(55) 



For large-amplitude driving, V3 ^Vc, one obtains 
i?3 ^ eVaFlGb + o 



(56) 
(57) 



Hence, as one would expect, large-amplitude driving is 
mostly sensitive to the behavior at large voltages. 

Another approach would be to DC bias the system to a 
given operating point Vq, and add a non-stationary term 
V3 cos(3a;oi) on top of this (this type of measurement 
is frequently used for the differential conductance). In 
this way, the measured third moment would be sensitive 
to that near the operating point only. In this case, the 
second moment would be of the form of Eq. H52|) with 
V{t) w Vq, only weakly dependent on V3 (except per- 
haps for Vo ~ Vc). On the contrary, the measured third 
moment would follow V'i, 

i?3 = e^F^GcV^ (58) 

where F3 is the Fano factor for the third moment and Go 
the differential conductance at the operating point. This 
approach would in principle work in any kind of nonlinear 
systems, provided that V3 is chosen appropriately. How- 
ever, this approach has the disadvantage that the second 
moment is relatively much larger than the third, making 
the averaging time long (see Sect.lIV|l. 



IV. LIMITATIONS OF MEASUREMENT 

In this section, we discuss the limitations in the mea- 
surement of higher moments of current fluctuations with 
the emphasis on the third moment. We consider the ef- 
fects of environmental noise added to the signal after am- 
plifying, the effect of a Gaussian noise in the electromag- 
netic environment, and the averaging time. 

A. Effect of amplifier noise 

In typical measurement systems, the signal coming 
from the filters has to be amplified before it is mixed and 
detected. This amplification introduces additional noise 
6Ia(t) in the output signal. In the following, we analyze 
the effect of such noise, assuming that it is Gaussian and 
uncorrelated with the fluctuations in the current coming 
from the sample. The second moment of such fluctua- 
tions is 

Saiio) = ^ [d{t~ t'){5Ia{t)5Ia{t'))e'^^'-''\ (59) 



whereas the average and the other odd moments vanish. 
Moreover, one has to assume that this noise is intrinsi- 
cally band-limited (due to finite i?C-times etc.), i.e., 

{5Ia{t)5h{t)) = -i= 1 duSaiuj) EE -^^"^ (60) 
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is finite. Due to this noise, the fluctuating signal before 
mixing is of the form 

X{t) = -^HBpit) * I(t) + SIa{t). (61) 

V 27r 



In the AC mixing scheme, the resulting signal after mix- 
ing is 

{Y) - {Y)o+-^ {htp * [fimBPit) * {im{6utr)] } , 

(62) 

where the first term has been calculated above and the 
second is due to excess noise. Such noise contributes only 
if the signal after band-pass filtering has a finite average 
and we apply f(t) — cos{LUQ{t + to)), 



(y) = (y)o + — ^tVicos(^oio). 



(63) 



For f{t) — cos{3u!Q{t + to)), the amplifier noise contribu- 
tion vanishes. 

In the AC statistics scheme, the amplifier noise does 
not contribute to the average signal. The second moment 
becomes 

{{5Yf) - {{dYf)o + -^Sa{iUo)iOm, (64) 

V 27r 



where the first term is that calculated above, and the 
second is due to amplifier noise within the band lu G 
[uq — ujjnjUjQ + LOjji]. As in the AC mixing scheme, the 
effect of amplifier noise is finite only for a non- vanishing 
average signal, 

{{SYf) = {{5Yf)Q + -^^a(^o)^™(/i) cos(a;oto). (65) 
47r 

Apart from the noise in the amplifiers between the fil- 
tering and mixing, there may be other sources of noise. 
The most important noise source, that in the electromag- 
netic environment seen by the sample is treated below. 
Another source of noise, possibly relevant in the mix- 
ing schemes, is the noise added by the mixers. However, 
there the dependence of the signal on the phase to can be 
used to distinguish the measured signal from the noise. 



B. Linear system with environment 

The effect of electromagnetic environment becomes the 
more important the higher moments one measures. For 
the second moment, the noise in the environment simply 
adds to that in the measured sample and scales the signal 
from the sample through the impedance ratio. Recently, 
Beenakker, Kindermann and Nazarov2^ showed that the 
environment has a strong effect on the third moment, as 
this includes the effect of current noise in the environment 
biasing the sample. In the following, we generalize their 
approach to the case of nonstationary driving. We also 
address the methods for filtering in these devices and 



is(t)Q !cjio(t) V+AV 



FIG. 3: Circuit studied for the effect of fluctuations in the 
electromagnetic environment on the third moment of the mea- 
sured fluctuations. Circuit consists of the measured system 
with impedance Zs and fluctuating current is(t), coupled to 
its environment with impedance Zo and Gaussian current fluc- 
tuations ioit). The arrow on the bottom of the figure indicates 
the chosen reference direction for the current. 



show how this yields the results of the general analysis 
in Sect.ini 

To properly deal with the filtering, let us assume that 
the measured system can be separated in two parts (see 
Fig.|2J): measured sample with impedance and a band- 
pass filter with impedance Zq{uj). The impedances of the 
measuring lines, current amplifiers, etc. can be included 
in Zq[lij). For simplicity, we assume that within the con- 
sidered frequency interval, Zg is frequency independent, 
and Zq{uj) is given by 



Zo{uj) = 



Ze, otherwise. 



(66) 



Here, Z^ is either much larger (in the case when one 
would measure directly the current fiuctuations) or much 
smaller (when converting the current fluctuations to volt- 
age fluctuations) than the sample impedance. For sim- 
plicity, all impedances are chosen real. Moreover, assume 
that the current fiuctuations in Zq (w) are Gaussian with 
the second moment given by 



Zo[uj) 



(67) 



Here Tq characterizes the noise temperature of the envi- 
ronment and in general it may differ from the tempera- 
ture of the sample. 

The intrinsic current fluctuations SI(uj) in the two 
impedances cause voltage fluctuations AV{oj) over them. 
As a result, the total current fluctuations at the two 
impedances are of the form 



A/o(t^) - SIo{lo) - 
A/,(w) = 6Isiuj) + 



AV{u;) 

Zo{iLl) 

AV{uj) 



(68) 
(69) 



Assuming low enough frequencies (compared to the fre- 
quencies describing charge relaxation), the two currents 
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should equal. This requirement leads to 
ZsZq{uj) 



AV{u;) = 
A/(c.; 



z. 



(70) 



SUuj) + M^SIo{^). (71) 
Z(uj) Z(uj) 



Here we defined Z{lu) = Zs+Zq{lj). The second moments 
of current fluctuations in the total system and voltage 
fluctuations over the impedances are given b}i2i 

Si{uj,uj') = {AI{u;)AI{lj')) 



Z{uj)Z{uj') 



Zo{u)Za{uj') f 

-bj[UJ,Uj) + —————bj[UJ, 



Z{lu)Z{lo') 



(72) 



5v(ijJ, uj') 



{AV{lu)AV{uj')) 
Z^Zoiuj)Zo{Lu' 



Z{uj)Z{uj') 



(73) 



where Sj{uj,uj') = {5Is{uj)6Is{uj')) measures the second 
moment of current fluctuations in the sample. 

One can easily find that if the fluctuations are mea- 
sured from the current, the optimal situation for noise 
measurement is Zp <^ Z^ and Z^ 3> Zs whereas for the 
measurement of voltage fluctuations, the maximum sig- 
nal is obtained with Zp k, Zg and Zi^ <^ Zs^ In practice, 
when the signal power has to be compared to the resolu- 
tion of the meters, the latter case is advantageous. 

The major invention by Beenakker et al. was the fact 
that as the magnitude of current fluctuations in meso- 
scopic samples depends on the amplitude of driving, the 



voltage fluctuations in the environment tune these fluc- 
tuations, and thereby induce cross-correlations between 
the current fluctuations in the sample and in the environ- 
ment. Following their work, we expand the second mo- 
ment of fluctuations, S1{V) in the driving voltage near 
the average voltage. 



5Is{u:)5Is{uo') = / dte'("+'^')*(<5/,(t))2 



dV 



(74) 



The first-order fluctuation term is enough to describe 
the environmental effect on the measured third moment. 
Now assume the driving voltage has the form given by 
Eq. © and expand the above two terms in Fourier har- 
monics. 



Sf{{V{t)))^Y.^r.cos{nu;ot) (75) 

ri=0 

dSf{V), ^ , , 

-\v={v{t)) = ^Dncos(nu;ot). (76) 



dV 



71=0 



For example, in a linear system in the limit eV 3> fc^T, 
driving with V{t) = VJj cos(3a;oi) yields 



Di = 0, Ds 



4^26 



(77) 



where Fo is the Fano factor for the second moment. 



Combining, we get 

6Is{uj)SIs{uj') ~ ^ X! \-'^n{S{uJ + oj' + jiujo) + S{uj + uj' - nujo)) + D,i{AV{uj + uj' + ncjo) + AV{uj + uj' - nujo))] ■ 

n 

(78) 

The third moment of the current fluctuations AI{uj) contains three parts: first coming from the "intrinsic" third 
moment of the measured sample, second from the cross-correlations between the sample and its environment, and 
the third from the current fluctuations of the sample driving the sample itselS^ 



Ri{u;,cu',ij") = {AI{lj)AI{uj')AI{lj")) = 
■Zo{uj) , Zoiu;') , Zo{uj") 



Z^ 



+ ksTo 



Z{uj) Z{oj') Z{uj") 



Z{uj)Z{lj')Z{lj") 



— ^ DnS„ 



n,a— ±1 

/ Zo{iLj + mPujo) Zo{uj' + mPuJo) Zo(w" + to/Jcjo) 
V Z{u; + m/3wo) Z{uj' + m/Suo) Z{uj" + m/Stoo) 



(79) 



d{ujs + {na + m(3)u}t)) 



Q,/3 = ±l 



where LOg = w -|- + lo" . This correlator is essential if one is able to directly measure the current fluctuations. In 
this case, specifying Z^^ Zg^ the fluctuations outside the band lo^lo\lo" £ ±[wo — Suj,ujo + Su] can be neglected and 
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within this band, 



{Z, + Zpf 



Ri iDi 



{2kBTo - ZsSo) - ^ DnSrn ) S{lJs ± c^o) 



2 {z, + Zpy 



m^O,n±m— 1 



(80) 



Other Fourier components of R and D vanish due to the filtering provided by Ze- 

Another important observable is the third moment of voltage fluctuations over the impedance Zq. This is given by 



P , , Z^Zo{lo)Zo{lo')Zo{^") 



Z{uj)Z{uj')Z{uj") 



R] + ZskeTo 



1 



1 



Z{uj) Z{uj') Z{lu") 



n,a— ±1 



Zo(ti' + to/Jwq) Zo{ll!' + miSuJo) Zq{lu" + mPujo) 
Z{uj + mPuo) Z(io' + mPuJo) Z{uj" + m/3wo) 



5(tJs + (na + m/3)a;o) 



a,/3=±l 



(81) 



Now the pass-band may be defined through Zp « Zs and the stop-band by Z^ <C Zs- Assume for simplicity that 
Zp = Zg. Then, if all the frequencies are within the pass-band, one obtains 



Rv{uJ,uj',uj") 



Z'i 



^5{ujs ± Wo) - -^'5(tJs ± 3wo) -I- -{ZsSo + 'iksTo) (DiS{ujs ± ujo) + D-iS{ujs ± 3wo)) 



r 



(82) 



In the stop-band, Ry is negligibleiSI 

As an example, consider a linear system with Fano 
factors F2 and ^3 (for the second and third moments) in 
an impedance-matched environment Zp = Zs and driven 
with V{t) = VsCosiSuJot), eVs > ksT. The third mo- 
ment of voltage fluctuations over the sample within the 
pass-band is 



Rv = Zt 
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—F,eh + -F,e^ 



(83) 



This equation illustrates how the measured third moment 
consists of the "real signal" (dependent on F3) and the 
parts coming from the sample noise driven by its own 
fluctuations (second term) and those in the environment 
(third term). 

In their paper^ Reulet et al. showed that the environ- 
mental effect can be accurately described by this theory 
after including the finite propagation time r in the coax- 
ial cable between the filter and the sample (see Fig. 
Their analysis is analogous in the case of slow driving, 
and indicates that a long cable can be used to decrease 
the effective noise temperature Tq of the environment. 



C. Averaging time 

In Ref. 20, the required time for signal averaging was 
several hours for each plot. Naturally, for any reason- 
able measurement, this time should not be much longer. 



To get an estimate of the required measurement time 
with the cyclostationary driving, we consider the AC 
statistics scheme (the requirements are analogous for 
the first scheme, only the prefactors of the expressions 
may slightly differ). There, one obtains a single output 
value in a time determined by the data integration time 
^mcas = Stt/wto. Assumc One has measured n values of 
current. The squared deviation between the measured 
value of the unbiased estimator ^3 for the third moment 
of the data and the true third moment can be estimated 
through the variance of k^i^ In the limit n 3> 1 it is given 
by 



1 



var(fc3) = -(9^2 - Ai3 - 6//2Ai4 + Me), 



(84) 



where /i^ are the z'th central moments of the true current 
distribution. Now let us require that the relative error is 
smaller than some percentage p of /^3, 



var(fc3) < p^^l\. 



This leads to the requirement 



1 



n > 



9^4 



1-6 



(85) 



(86) 



From this example, it is thus clear that any additional 
noise from the measuring setup increases the measure- 
ment time through the increase of the even moments. 

In the typical limit where the current / through the 
sample exceeds e times the bandwidth Ud (e.g., for / = 
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100 nA, J/e = 620 GHz), it is safe to neglect the cumu- 
lants of higher order than three, i.e., the measured data 
is almost Gaussian. In this case the above requirement 
reduces to 

6 ^il _ 327r {2So + S2cos{2ujoto)f 



n > 



Sp^ujd (3i?i cos{ujoto) 



i?3 cos(3t^oio))^' 
(87) 

where the last form was obtained using Eqs. I|29l3n|l . 

To get an idea of the required measuring times, con- 
sider a linear system driven with the voltage V3 cos(3cL;oi) 
in the shot-noise limit eVs 3> ksT, and in the case when 
the sample noise dominates the second and third mo- 
ment of the fluctuations (see Subs. UlL^ . In this case 
the required number of points for = is 



n > 



2048 Fi I3 



Using the values I3 = 100 nA, = ^^d = 27r(10MIIz), 
and p = 0.05 for a tunnel junction (F2 — ~ 1), we 
obtain 71 > 2.7 • 10^ and T = 27rn/u;„ > 27 s. 

If the current is lowered, also the measuring time is 
lowered as the "skewness" A'.3//i2^^ increases with de- 
creasing average. Of course, the estimate is valid only 
in the shot-noise limit eV 3> fc^T and when the sample 
noise dominates the amplifier noise. Therefore, optimal 
signal is expected for the minimum values of current with 
which the signal moments are still determined by shot 
noise. Note that typically increasing the impedance of 
the measured sample decreases the possible bandwidth 
for impedance matching. However, also the shot-noise 
limit eV 3> ksT is obtained with a lower current. As a 
result, the overall averaging time for the third moment 
does not change much. 

If the amplifier noise 5j = 2fcsTo/Zo dominates the 
shot noise, we get in the otherwise same case as above 
the requirement 

20487T fkBToV h 
3p^ \eI3Z0J eujd 

Now, as an example, for Tq = I3 = l/iv4, Zq = 
50il, bJm = i^d — 27r(100MHz) and p = 0.05 the required 
number of measurements for a tunnel junction would be 
n > 4.4 • 10^° and T > 440 s. 

A similar estimate for the second moment would yield 
a sample size (in the limit n ^ 1) 

^ r(A*4-/i^) (90) 



n > 



and for the example considered above. 



1 fSn ^ 



0.17 

p2 



(91) 



independent of the current or the bands (in the Gaus- 
sian limit where the prefactor is slightly different, this 
is analogous to the Dicke radiometer formulai^Si) . Hence, 
for p = 0.01, we would need a few thousand samples, and 
the required measuring time with MHz sampling would 
be of the order of milliseconds. 



V. SUMMARY 

In this paper, we detail methods to measure the mo- 
ments of current fluctuations by employing a nonstation- 
ary driving signal. Such driving overcomes the require- 
ment for the wide-band measurements and gives a way to 
confirm the measurement results through the variation of 
the phase between the signal and the mixing. However, 
cyclostationary driving makes rise to higher harmonics 
of the signal, which have to be taken care of through a 
proper Fourier analysis. In Sect. IIIII we discuss a few 
examples of such an analysis. 

Section IIVI shows that the measurement of the third 
moment always involves the effect caused by the mea- 
suring setup. The averaging times discussed at the end 
of the same section indicate that this type of "classi- 
cal" measurements are practical for the measurement of 
the third and perhaps still for the fourth cumulant, but 
higher cumulants seem to be out of the time constants 
set for any reasonable project. Therefore, schemes where 
the detector is a mesoscopic system and placed near the 
sampl oi§ii3S i, overcoming the low-bandwidth restrictions, 
are preferable. 

Let us outline a practical scheme for the detection of 
the third moment/cumulant of current fluctuations with 
cyclostationary driving. Assume one measures a sample 
with (possibly voltage-dependent) impedance Zs- 

1. Construct the matching circuit that approximately 
matches the input impedance of the amplifier to 
that of the sample, and filters the signal from the 
narrow band u G ±[ujo — Suj, ujq+Suj]. Note that the 
measurement outcome will depend on the details of 
the matching, as discussed in Subs. HVBl 



2. Drive the sample with voltage V{t) — Vq + 
V3 cos{3ujot) where either Vq = or Vq ^ V3 ("dif- 
ferential measurement", see Subs. IhTHI . 



3. Either mix the signal twice with itself and once with 
f{t) = cos{3uJo{t + to)) fSubs. lnX2|) or only with 
f {t) = cos{uja{t + to)) fSubs. IIIB} . Average the 
signal (mixing scheme) or its third moment (statis- 
tics scheme) for time exceeding that indicated in 
Subs. IIVCI The outcome is given by Eq. or 
Eq. H30|) . where i?3 is taken from the cos(3(jJoO" 
component of the Fourier series of R{V{t)) as in 
Subs.imSl 
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APPENDIX A: MOMENTS IN THE AC 
STATISTICS SCHEME 



statistics scheme can be written as 



The contribution of the Pih harmonic of the n'th mo- 
ment, M^, to the n'th measured moment in the AC 



{iSYr)i = J dh... dt^H'^pih) . . . i7™p(i„)/(ti) . . . /(t„) j dt'HBPih ~t')... HBp{tn - t') cosM). (Al) 



r 



The convolution between the iJsp-functions filters out 
most of the harmonics I. It only leaves I = 1, 3, . . . , n in 
the case when n is odd and Z = 0, 2, . . . , n when n is even. 



Now, use the fact that HBp{t) = 2cos{ujQt)Hf^p{t), 
where -ff|,p is a low-pass filter with band lo e [—5uj, 5uj\. 
One obtains 



J 



{{5Yr)i = Ml 



^ j dh... dt^dt'H^pih) . . . H'^p{t^)Hip{h ~t')... Hlp{t^ - t')G{h, . . . , i„, t'), (A2) 



where 



G{ti, . . . , i„, t') = f{ti) . . . f{tn) cos(wo(ii - t')) . . . cos(wo(t„ - t')) cos{lujQt'). 

I 



(A3) 



If LUrmSuj < LOo/n, only the time-independent part of 
G{ti^ . . . ,tn, t') survives the remaining filtering. This can 
be computed from 



n+l /•2ir/tJo 



dti... dtndt'G(ti 



t t') 
(A4) 

However, G^^q may still depend on the relative phase 
to between the driving and the mixed function f{t) = 
cos(ajo(t + to))- Equivalently to above, G^q is only finite 
if I is one of the values that survives band-pass filtering. 
Moreover, it is straightforward to show that G'^'^(to) — 



gl^ cos{lu)[)tQ) /2^'^^ . The number g\-^ for some of the lowest 
(nontrivial) values of n, I is tabulated in Tabled 



n 


1 


2 


2 


3 


3 


4 


4 


4 


5 


5 


5 


6 


6 


6 


6 


I 


1 





2 


1 


3 





2 


4 


1 


3 


5 





2 


4 


6 


5n 


1 


1 


1 
2 


3 

4 


1 

4 


3 
4 


1 
2 


1 

8 


5 
8 


5 

Ifi 


1 

ifi 


5 
8 


15 

32 


3 

Ifi 


1 

32 



TABLE I: Values of for diflGrcnt moments n and different 
harmonics I. 

The remaining part of the integral in Eq. (|A2() is now 
a combination of low-pass filters: 



dhdh . . . dt,,dt'Hrp{h) . . . Hrp{tn)Hlp{t, -t')... Hlp{t,, - t'). (A5) 



Singling out one of the U integrals, it is straightforward 
to show that 



dUHl'p{U)Hip{U - t') = ^T^^ \H-"j:p{uj)Hip{Lo) 
= V2^Htp{t'), 

(A6) 



where H'[p{t') is a low-pass filter with bandwidth uj^ 
rmii{uj„i, Suj) {uj G [—aJd^LOj]). Thus we get 



= ^ (^) ^ Mlglcos{lu;oto)nnM (A7) 
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Here 



I 



(27r)"/2-i 

doJi... duJn~-lHip{uJi) . . . Hip{uJn-l)Hip{uJi + . . . UJn-l) 



n-1 

(27r)"/2-i- 



and a„ is a number whose values for lowest n are 



(A8) 



Finally the contribution of the Z'th harmonic of the n'th 
moment on the 71 'th moment of the measured quantity is 
given by 



= ^^r^5>nc^r cos(Zc^oio). (AIO) 



ai — 1 
"3 = 3 
as = 115/12 

= 17407/360 
ag = 18063361/40320 - 1440. 



a2 = 2 
a4 = 16/3 

ag = 88/5 (A9) The lowest six moments are listed in Eqs. (|28I33|) . 

ag = 53752/315 
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